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ABSTRACT 

The dwarf spheroidal galaxies (dSph) of the Milky Way are among the most attractive targets 
for indirect searches of dark matter. In this work, we reconstruct the dark matter annihilation (J- 
factor) and decay profiles for the newly discovered dSph Reticulum 11. Using an optimized spherical 
Jeans analysis of kinematic data obtained from the Michigan/Magellan Fiber System (M2FS), we 
find Reticulum IPs J-factor to be among the largest of any Milky Way dSph. We have checked 
the robustness of this result against several ingredients of the analysis. Unless it suffers from tidal 
disruption or significant inflation of its velocity dispersion from binary stars. Reticulum H may provide 
a unique window on dark matter particle properties. 

Subject headings: galaxies: dwarf — galaxies: individual (Reticulum H) — dark matter — gamma 
rays: galaxies — methods: statistical — stars: kinematics and dynamics 


1. INTRODUCTION 


Along with the Galactic center and galaxy clusters, 
the dwarf spheroidal galaxies (dSph) of the Milky Way 
have been identified as promising targets for indirect 
dark matt er (DM ) searches (see e.g. |Strigari||2013 Con¬ 


rad et al. 2015). Their low astrophysical background, 


high massAo-lig ht ratio, and proximity make th em com¬ 


pelling targets ( Lake| 1990 Evans et al. 2004). About 


twenty-five Galactic dSphs were known as of early 2015, 
and their observation by y-ray telescopes has thus far 
shown no significant emission, leading to stringent con¬ 
straints on 


tion cross-section dAcciari et al.|20I0 Paiano et al.|2011 

Abramowski et ai.||2l)14t 

Geringer-Hameth et ai.||2l)14 

h'ermi-LAi' Collaboration 

m5F^ 


has led to the discovery of nin e new potential Milky-Way 
satellites in the Southern sky (Koposov et al.|[2015 


DES 


Collaboration et al. 2015). The nearest object, Reticu- 
30 kpc), is particularly intriguing, as 
evidence of y-ray emission has been detected i n its direc¬ 
tion using the public Fermi-LAT Pass 7 data. [GeringerT] 
Sameth et al. (2015b) determined the probability of back¬ 
ground processes producing the observed Ret 11 gamma- 
ray signal to be between p = 0.01% and p = 1%, de¬ 
pending on the background modelling. An analysis of 
the new objects published simu ltaneously by |The Fermi-| 
LAT Gollaboration et al. (2015), based on the unreleased 
Pass 3 data set, reported no significant detecti on, though 
the strongest hint was for Ret 11 with p = 6%. Hooper fc 
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Linden (2015) subsequently performed a similar analysis 
with public Pass 7 data, finding a p value of 0.16%. 

In any case, a robust determination of Ret IPs DM con¬ 
tent is crucial in order to constrain particle nature of DM. 
Reticulum H was found to be a DM-dominated dSph 


of 

Walker et 

al. 1 

20151) 

posov et al. 1 

2Ul5l 

M. H 


Simon et al. (2015) and 


Ko- 


ere, we reconstruct the DM an- 


a spherical Jeans analysis applied to stellar kinematic 
data obt ained with the Mich igan/Magellan Fiber System 
(M2FS) ( Walker et al.|201^. We use the optim ized Jeans 


analysis setup from Bonnivard et al. 
pute the astrophysical J- and D-tac 


( |2015a|b[ ), and com- 
ors, for annihilating 


and decaying DM respectively, from the reconstructed 
DM density profiles. We cross-check our results by vary¬ 
ing different ingredients of the analysis and evaluate the 
ranking of Ret H among the most promising dSphs for 
DM indirect detection. 

2. ASTROPHYSICAL FACTORS, JEANS ANALYSIS AND 
DATA SETS 

2 . 1 . Astrophysical factors 

The differential y-ray flux coming from DM annihila¬ 
tion (resp. decay) in a dSph galaxy is proport ional to the 

so-called ’astrophysical factor’ J (resp. D) (Bergstrom 

- 


J= did (resp. D = 


=JJpDM{l,^)dld^^, ( 1 ) 


which corresponds to the integration along the line-of- 
sight (l.o.s.) of the DM density squared (resp. DM den¬ 
sity) and over the solid angle AD = 27r x [1 — cos(aint)], 
with ttint the integration angle. This quantity depends 
on both the extent of the DM halo and the mass density 
distribution, and is essential for constraining the DM par¬ 
ticle properties. All calculatio ns of astrophysical factor s 
are done with the CLUMPY code (Charbonnier et al.|2012 ), 
a new module of which has been speciticaily developed 
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Bonnivard et al. 


to perform the Jeans analysifQ 

2 . 2 . Jeans analysis 

Several approaches have been developed to infer the 
DM density 


Battaglia et al. 

2013 

Strigari 

2013 

Oswalt 


of astrophysical factors dStrigari et al.||2007 Essig et al. 

20101 ICharbonnier et al.| 2011 

Cholis fc 

Saluccl 2012 

Ceringer-Sameth et al. 2015al 

Bonnivard 

"eTaT 2015b|). 

We reter the reader to Donnivard et al. 

(|20i5a|) tor a 


work. Here, we summarize the main ingredients. 

Assuming steady-state, spherical symmetry, and neg¬ 
ligible rotational support, the second-order Jeans equa¬ 
tion, o btained from the collisionle ss Boltzmann equation, 
reads (Binney fc Tremaine||2008 1: 


^ d , 2\ I 

V dr 


{r)vl 


GM{r) 


( 2 ) 


with u[r) the stellar number density, u^(r) the radial 
velocity dispersion, j5a.m(r) = 1 — the velocity 

anisotropy, and M{r) the mas^ enclosed within radius 
r. After solving Eq. (§ and projecting along the l.o.s., 
the (squared) velocity aispersion at the projected radius 
R reads 


^liR) = 


m) 


l-/3ani(T’)- 


^ i/(r) v^{r) r 


dr, (3) 


with E(i?) the surface brightness profile. We compare 
the l.o.s velocities of the stars to the projected veloc¬ 
ity dispersion Up, computed using parametric forms for 
the unknown velocity anisotropy /3ani('c) and DM density 
profile Pniufr). We use the following likelihood function 
( Strigari et al.|[2007 ) 






(27r) 


- 1/2 


=i^cr2(i?.)_|_A2_ 


■ exp 


{Vj-vY \ 

al{Ri)+IS.lJ 


(4) 


which assumes a Gaussian distribution of l.o.s. stellar 
velocities Vi, centered on the mean stellar velocity v, 
with a dispersion of velocities (at the radius Ri) com¬ 
ing from both the intrinsic dispersion ap{Ri) and the 
measurement uncertainty Ay.. Probability density func¬ 
tions (PDFs) of the anisotropy and DM parameters are 
obtained with a Markov Chain Monte Carlo (MCMC) 
engin^ and are used to compute the median and cred¬ 
ible intervals (CIs) of the astrophysical factors for any 
integration angle. 


in 

Bonnivard et al. 1 

2015a), the DM density is 

by 

an Finasto pro 

tile (Merritt et al. 

2006 

I, 


aniso tropy and light profile s are given by Baes fc van 
Hese ( Baes fc van Hese[2 007 ) and Zhao-Hernquist (Hern- 


quist 


1991) Zhao ' 1990) parametrisations, respectively. 


° This upgrade will be publicly available in the new version of 
the software (Bonnivard et al., in prep.). 

^ The mass is dominated by DM, and we neglect the stellar 
component. 

^ We use the GreAT toolkit ( Putze|2011 Putze &: Derome|2014 1. 
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Figure 1. Projected stellar density profile of Ret I I. derived from 
the photometric catalog of |Koposov et al.| l|2015ll. Overplotted 
(red line) is the best-fitting model (we note mat the fit is to the 
unbinned data), which is the sum of cont ribu tions from Ret II itself 
and a constant background (see Section |2.3| |. Dotted lines enclose 
68% CIs for the projection of !^(r). 


The large freedom allowed by these parametrisations was 
found to mitigate possib le biases of the Jeans analysis 
( Bonnivard et al.]|2015a ). Finally, the extent of the DM 
halo is computec T using t he tidal radius estimation as in 
Bonnivard et al. (2015b). 


2.3. Data set 

Surface brightness data — We fit the stellar number den- 
sity profile i/jr) o f Ret I I following the procedure that 
Bonnivard et al. (2015b) use for ‘ultrafaint’ dSphs (see 
their section 3.1). We consider a flexible Zhao-Hernquist 
model for the 3D profile. 


,Zhao 


(r) = 


(r/r*)T'[l -h (r/rj)"]!/^-^)/" 


(5) 


where the five parameters are the normalization v*, the 
scale radius r*, the inner power law index 7 , the outer 
index /3, and the transition parameter a. Along with 
an additional free parameter Ebkd that represents a uni¬ 
form background density, these parameters then specify 
a model for the projected stellar density: 


r 

Smodel(J?) = 2 / 
JR 


/{r)7 


\/r^ — 


dr + Ebkd- 


( 6 ) 


W e ht this model to the photometric catalog generated 
by Koposov et al. (2015), which provides positions, col- 
ors, and magnitudes of individual stars detected as point 
sources. From the raw catalog, we first identify pos¬ 
sible members of Ret H as point sources (selected as 
sources with Sextractor ‘spread’ parameter < 0.01 in 
the g-band) whose extinction-corrected g — r colo rs place 
them within 0.25 dex of the Dartmouth isochrone (Dotter 
et al.||2008|), calculated for a stellar population with age 
12 Cyr, metallic ity [Fe/H] = —2.5, an d distance modulus 
m — M = 17.4 (Koposov et al. 2015). To the unbinned 
distribution of projected positions tor the N = 12470 
RGB candidates identified within 1.5° of Ret H’s center, 
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Figure 2. Top: velocity dispersion profile of Ret II and recon¬ 
structed median and credible intervals (solid and dashed black lines 
respectively), as well as best fit^ (long dashed red lines). Bottom: 
distribution of membership probabilities as a function of the pro¬ 
jected radius R and the departure from the mean velocity (z-axis, 
blue to red color) for the eighteen stars with 0. The size of 

the points is proportional to the velocity uncertainty. See text for 
discussion. 


we fit 2D projections of v{r) according to the likelihood 
function: 

N 

^2 OC R I]i„odel(Ri)- (7) 

i=l 


As in Bonnivard et al. (2015b ), the fit is done with the 
software package IVIultiNest (Fbroz fc Hobson 2008 Feroz 
et al. 2009 2013), and we use the samples from the poste- 
rior ■PUFV to propagate the IDht profile uncertainty into 
the Jeans analysis. Figure shows the ht to the pro¬ 
jected stellar density profile of Ret II (dashed red line), 
with the contributions from Ret II itself and from the 
constant background (solid black and blue lines respec¬ 
tively). 

Kinematic data — We use the Ret II stellar kinematic 


data set from Walker et al. (20151, obtained with M2FS. 
It consists of projected positions and l.o.s. velocities for 
38 individual stars, as well as an estimation of their mem¬ 
bership probability Pi. The latter, obtained using an ex - 
pectation maximization algorithm (Walker et al. 2009), 
quantifies the probability that a given star belongs to the 
dSph or to the Milky Way foreground. 

The top panel of Figurel^ presents the velocity disper¬ 
sion profile of Ret II, as well as its reconstruction with the 
Jeans analyst^ The bottom panel of Figure]^ shows the 
distribution ofmembership probabilities as a function of 
the projected radius R and the departure from the mean 
velocity (color-coded), for stars with non-zero Pi. As 


Table 1 

Astrophysical factors for Ret II {d = 30 kpc). For five different 
integration angles, the median J (resp D)-factors as well as their 
68% and 95% CIs are given. Note that possible triaxiality of the 
dSp h galaxies adds a system atic uncertainty of ±0.4 (resp. ±0.3) 
jBonnivard et al.|2015all and is not included in the quoted 
intervals. 


aint 

loglo(.I(«int)) 

loglo(t^(«int)) 

[deg] 

[ J/GeV^ cm-5p] 

[D/GeVcm-2p] 

0.01 


.. 7-h0.6(-hl.oH 

^'^•'-0.3(-0.5) 

0.05 

.0 o+0.5(+l.l) 

7 -h0.5(-hl.0) 
•^_0.3(_o.6) 

0.1 

1« o+0.6(+1.2) 

-Lo-0_o.5(-0.8) 

-,7^-h0.6{-hl.l) 

•'^-0.4(-0.6) 

0.5 

+ 1.0(+1.7) 

-Ly-o_o.7(-1.3) 

-h0.7{-hl.2) 

18.8_0 7(_i 1 ) 

1 

19 

1 Q q+0-8{-|-1.4) 

J-y-'5 0.9(-1.4) 


="1 GeV^ cm-5 = 2.25 X IQ-’^ kpc”® 
*^1 GeVcm“^ = 8.55 x 1O~^®M0 kpc“^ 


pointed out in Bonnivard et al. (2015b), a large fraction 
of stars with both intermediate Pi (U.l < Pi < 0.95) and 
large departure from the mean velocity hints at Milky 
Way foreground contamination, which can affect the J- 
and D-factor reconstruction. For Ret II, only one star 
shows an interme diate Pi (Ret2-142 in the catalog of 
Walker et al.|2015 with Pi = 0.69), with a very small de- 
parture from the mean velocity. Therefore we do not ex¬ 
pect a strong sensitivity to fore ground contamination. I n 
this study, and as advocated in Bonnivard et al. (2015b), 
we use the data with Pi > 0.95 (sixteen likel y members, 
one less than identihed by |Walker et al.|2015 after exclu¬ 
sion of Ret2-142) as our fiducial setup. 

3. RESULTS 

Figure displays the J- (top) and D-factors (bottom) 
of Ret Unreconstructed from the Jeans/MCMC analy¬ 
sis, as a function of the integration angle Oint- Solid lines 
represent the median values, while dashed and dash-dot 
lines symbolize the 68% and 95% CIs respectively. Our 
data-driven Jeans analysis gives large statistical uncer¬ 
tainties due to the small size of the kinematic sample, 
co mparable to those obtaine d for other ‘ultrafaint’ dSphs 
by Bonnivard et al. (2015b) (see also Figure |4|). Table 
[B summarizes our results for the astrophysical factors of 
Ret II. 

We cross-check our findings by varying different ingre¬ 
dients of the Jeans analysis. The resulting J-factors are 
shown in Fi gure First, we perfor m a binned Jeans 


2015b) of the kinematic 
ts. Second, we cal culate 


analysis (see Bonnivard et al. 
data, and find compatible resii 
the b ootstrap mean and dispersion of the J-factor (|Efron 
|1982| ). For this purpose, we generate 500 bootstrap re- 
samplej^ by drawing with replacement 16 stars among 
the 16 of the original sample with Pi > 0.95. The re¬ 
sults are in excellent agreement with the MCMC analy¬ 
sis. Finally, we use all 38 stars of the sample but weight 
the log-likelihood fu nction of Eg. (1^ by the member¬ 
ship probabilities Pi (Bonnivard et m. 2015b). As only 
one star shows an intermediate membership probability 


® The binned data and associated velocity dispersion reconstruc¬ 
tion are only shown for illustration purposes. The final results are 
obtained with an unbinned analysis. 


The best-fit DM profile and anisotropy parameters for each 
sample are obtained by maximizing the likelihood of Eq. J- 

factors were then computed for these best fitting profiles. 
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Figure 3. Median (solid), 68 % (dashed), and 95% (dash-dot) CIs 
of the J- (top) and D-factors (bottom) of Ret II, as a function of 
integration angle, reconstructed from our Jeans/MCMC analysis. 


0.1 < Pi < 0.95, we obtain very similar results. These 
two tests confirm that the reconstruction of the astro- 
physical factors of Ret II is not significantly affected by 
outliers. This is not always the case, notably for Segue I 
(Bonnivard, Ma urin fc Walke r , in pr ep.). 

We note that Simon et al. (20151 independently per¬ 
formed an analysis of the iVI2h'S Ret II spectroscopic data 
and found a slightly smaller J-factor. This can be traced 
to their choice of priors and light profile (L. Strigari, 
private communication). A detailed comparison will be 
presented in Geringer-Sameth at al. (in prep.). 
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Figure 4. Comparison of the J-factors at aint = 0.5° obtained 
for R et II (red circle) a nd for the potentially brightest objects 
from |Bonnivard et al.| ||2015b|| (blue squares), with the same 
Jeans/MUMU analysis. Ket 11 is comparable to Wil I in terms 
of J-factors, but slightly below Coma and UMa IL A 0.4 dex sys¬ 
tematic uncertainty was added in quadrature t o the 68% CIs to 
accoun t for possible triaxiality of the DM halo ||Bonnivard et aL| 
|2015a||. Also shown are the J-factors obtained forget 11 by varying 
ditlerent ingredients of the analysis - see Section]^ 


certainties in the J-factors. Moreover, the Jeans analysis 
assumes all of these objects to be in dynamical equilib¬ 
rium, but tidal interactions with the Milky Way could 
artificially inflate the velocity dispersion and therefore 
the astrophysical factors. UMa II, and to a lesser ext ent 


Coma, appear to be experiencing tidal disturbance ( Si- 
mon Sz Gehal 120071 IFellhauer et al.||2007l | Munoz et al. 


Mop [Smith et al. 


IFellhauer et al. [2007[ [M- 
20I3D, while Wil I may 


may show non¬ 


equilibrium kinematics (Willman et al. 2011). Caution 
is therefore always advised when interpreting the astro- 
physical factors of these objects. The dynamical status 
of Ret II is not yet clear. Its flattened morphology may 
signal ongoing tidal disruption. However, the available 
kinematic data do not exhibit a significant velocity gra¬ 
dient that might be assoc iated with tidal streaming mo¬ 
tions ( Walker et al.||2015 |. 


4. COMPARISON TO OTHER DSPHS 


one other dSphs in 

Bonnivard et al. ( 

2015b 1. In Figure 

3 

we compare the J-factors (tor aint = 

0.5") ot Ret H to t 

he 

brightest objects identified in 

Bonnivard et al. 

(2015b[^^ 


J-factor, but slightly below Coma Berenices and Ursa 
Major II. Its CIs are typical of an ‘ultrafaint’ dSph, and 
significantly larger than the uncertainties of ‘classical’ 
dSphs. 

Interpreting the possib le 7 -ray signal in Ret II in 
terms of DM annihilation (Geringer-Sameth et al.|2015b 


Hooper fc Linden||2015 ), one would expect similar emis¬ 


sions from the dSphs with comparable J-factors, such as 
UMa H, Coma, and Wil I. However, no excess was re¬ 


ported from these latter objects ([Gerin ger-Sameth et al. 
2014 Fermi-LAT Collaboration |2(JI5| ). This could be 
explained by the large statistical and systematiej^ un- 


5. CONCLUSION 

We have applied a spherical Jeans analysis to the newly 
discovered dSph Ret H, using s ixteen likely member s 


from the kinematic data set of Walker et al. (|M^. 
We einpl oyed the optimized setup of [Bonnivard et al. 
(2015a|b), which was found to mitigate several biases 


of the analysis, and checked that our results are robust 
against several of its ingredients. We find that Ret H 
presents one of the largest annihilation J-factors among 
the Milky Way’s dSphs, possibly making it one of the 
best targets to constrain DM particle properties. How¬ 
ever, it is important to obtain follow-up photometric and 
spectroscopic data in order to test the assumptions of dy¬ 
namical equilibrium as well as to constrain the fraction of 
binary stars in the kinematic sample. Nevertheless, the 
proximity of Ret H and its apparently large dark matter 
content place it among the most attractive targets for 
dark matter particle searches. 


1^1^ Segue I may have a highly uncertain J-factor (Bonnivard, 
Maurin & Walker, in prep.). We show it only for illustration pur¬ 
poses. 

1^2 The latter comes from a possible triaxiality of th e dSph (0.4 
and 0.3 dex for annihilation and decay respectively, see|Bonnivard| 
[et al.|2015a[ |, and depends on the l.o.s. orientation with respect to 
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